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Abstract 

We introduce several bilocal algorithms for lattice self-avoiding walks that provide 
reasonable models for the physical kinetics of polymers in the absence of hydrody- 
namic effects. We discuss their ergodicity in different confined geometries, for instance 
in strips and in slabs. A short discussion of the dynamical properties in the absence 
of interactions is given. 



1 Introduction 



The lattice self-avoiding walk (SAW) is a well-known model for the critical behaviour of 
a homopolymer in a solvent and it has been extensively used in the study of various 
properties of heteropolymers 0,^. Experiments are usually performed using monodisperse 
solutions and thus extensive work has been done to devise Monte Carlo algorithms to 
simulate fixed-length SAWs. Historically, the earliest algorithms used a local dynamics 
0: at each step a small part of the walk (usually 2-4 consecutive beads) was modified. 
Although easy to implement, these algorithms suffer a very serious drawback: as shown 
by Madras and Sokal 0], any local algorithm is not ergodic and simulations span only 
an exponentially small part of the phase space. A different algorithm was inspired by an 
attempt to model the true dynamics of the polymer in the solvent: the reptation algorithm 
|7|-pi]|. However, it was soon realized l^,^ that it is not ergodic because of the possibility 



of configurations with trapped endpoints. These ergodicity problems can be solved using 
chain-growth algorithms [p!IH l3[| or non-local algorithms |1I^-[18|. In the absence of any 
interaction, non-local algorithms are very efficient. For instance, in the pivot algorithm 
the autocorrelation time for global observables increases linearly with the number of steps 
N, which is the optimal behaviour since it takes a time of order simply to write down 
the walk. 

Non-local algorithms are extremely efficient in the absence of interactions and in free 
space. However, in the presence of strong attraction and in finite geometries non-local 
moves are largely inefficient^. Indeed, if the interactions are strongly attractive (for instance 
at the 6'-point or in the collapsed phase), typical configurations are compact so that the 
probability of success of non-local moves is very small. Therefore, the algorithm becomes 
inefficient. In the presence of surfaces, non-local algorithms are not even ergodic in general. 
For instance, as we will show, the pivot algorithm is not ergodic in a strip. Moreover, even 
when they are ergodic, they are not suited to study surface transitions since non-local 
moves will generate new walks with large energy differences and thus they will be rejected. 
Again the dynamics will be very slow. 

Stochastic algorithms are also important for the study of the physical kinetics of poly- 
mers. There is now a widespread interest in understanding the dynamics of homopolymers 
in the collapsed phase |^ and of heteropolymers near the folding temperature p2|-p^. In 
these studies one should use a dynamics which resembles the true dynamics of the molecule 
in a solvent. Local algorithms cannot provide good physical models because, as we said 



^An exception is the class of moves introduced in [1^ that do not change the position of the beads but 
only the connectivity of the walk. However these moves do not change global size observables and are of 
interest only for the study of maximally compact configurations. A general lower bound on the efficiency 
of non-local algorithms with a Metropolis test is given in Ref. EG] . 
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above, they are not ergodic and because, as noted in Ref. [^, they do not allow diffu- 
sive motions along the chain: for instance, they cannot move around assembled pieces of 
secondary structure^... As we said above, the ergodicity problem can be solved by us- 
ing non-local algorithms, and indeed, non-local moves have been considered in p5|, p3| , [3^ . 
However a non-local dynamics which involves rigid deformations of a large section of the 
polymer is unphysical, and therefore cannot give realistic results for the physical kinetics. 

In this paper we wish to discuss a family of algorithms that use bilocal moves: a bilocal 
move alters at the same time two disjoint small groups of consecutive sites of the walk 
that may be very far away. Since a small number of beads is changed at each step, these 
algorithms should be reasonably efficient in the presence of interactions, and thus they 
can be used in the study of the collapsed phase and of the folding of heteropolymers. 
They generalize the reptation algorithm and use a more general class of moves that was 



introduced by Reiter [^. Similar moves were introduced in Ref. ||36| and were applied to 
the study of ring polymers^. 

We will study in detail the ergodicity of these algorithms and we will show that, with 
a proper choice of moves, they are ergodic even in some constrained geometries, e.g. in 
strips, slabs, and generalizations thereof. These results have been obtained for SAWs with 
nearest-neighbour jumps on a (hyper-)cubic lattice. However, they can be easily generalized 



to different lattice models, for instance to the protein model proposed in Ref. Since 
these algorithms change a small number of beads at each step, they can be used to study the 
physical dynamics of dilute polymers in the absence of hydrodynamic interactions and are 
therefore a suitable generalization of the Rouse model [^,0 in the presence of excluded- 
volume effects. Therefore they provide a good model for the collapse of homopolymers and 



^In spite of these problems a local dynamics is of widespread use, see e.g. [p6|-p2[. Ref. claims that, 
at the presently investigated values of N, non-ergodicity effects should be small. A quantitative study of the 
Verdier-Stockmayer algorithm in two dimensions Q was presented in Ref. |6j . The percentage of walks that 
do not belong to the ergodicity class of the straight rod is indeed small, precisely 0.0067%, 0.0061%, 0.0041% 
for N = 11, 13, 15. However, it should be noted that this percentage increases with N (each ergodicity 
class contains only an exponentially small fraction of the walks ) and that the walks that do not belong to 
the ergodicity class of the straight rod correspond to compact configurations: therefore, larger systematic 
deviations are expected in the collapsed regime. If we indicate with CN{n) the number of walks with N steps 
and n nearest- neighbor contacts and with dN{n) the corresponding number of walks that do not belong 
to the ergodicity class of the straight rod, a good indication of the deviations expected in the presence of 
strong attraction is given by Rn = dN{nmax)/cN{nmax) and 5Ar(e) = J2n^"'''^Nin)/ J2n^"^'^Nin)- Here 
"iT-max is the maximum number of possible contacts for a given N] clearly Rn — Sn{oo). For the Verdier- 
Stockmayer algorithm we considered above, we have Rn — 3.2%, 1.4%, 5.8% for N = 11, 13, 15, and 
Sn{^) — 0.20%, 0.20%, 0.16% for the same values of N. Clearly, the systematic error is not completely 
negligible in the compact regime. 

•^However, it should be noted that the algorithm of Ref. [Bm is not ergodic. 
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for the folding dynamics^. 

The paper is organized as follows. In Sect. 2 we introduce several local and bilocal 
moves and we define three bilocal algorithms. In Sect. 3 we discuss their ergodicity, 
determining the minimal set of moves that make each algorithm ergodic. This is important 
in order to understand the dynamical behaviour. Indeed, algorithms that are ergodic only 
if rarely accepted or rarely proposed moves are included, are expected a priori to have a 
slow dynamics. Sect. 4 contains a detailed presentation of their implementation. In the 
last Section we present a brief discussion of the expected dynamic critical behaviour in the 



absence of interactions. A detailed numerical study will appear elsewhere |40 



2 Definition of the algorithms 

In this paper we will consider SAWs with fixed number of steps and free endpoints in 
finite geometries. 

More precisely, we consider a (i-dimensional hyper-cubic lattice and define the following 
set of lattice points: given an integer D such that 1 < D < d — 1, and {d — D) positive 
integers wd+i, . . . , Wd, we define Cd{wd+i, • • • , Wd) as the set of lattice points (ni, . . . , n^). 
Hi G Z, such that < rij < Wi, foi i = D + 1, . . . ,d. We will callQ Coiwo+i, • • • , Wd) a 
L)-dimensional cylinder. If D = {d — 1), we will speak of a strip if d = 2 and of a slab if 
d = 3. The number Wi will be called the width of the cylinder in the i-th direction. Note 
that we will always assume D > 1, so that at least the first direction is infinite. 

We will then consider SAWs of length A^ confined inside a cylinder. A SAW a; is a set 
of A^ + 1 lattice points uj^O), . . ., uj{N) such that: uj{i) and u^i + 1) are lattice nearest 
neighbours; uj{i) ^ u{j) for any i ^ j; u{i) G Cd{u!d+i, ■ ■ ■ ,Wd) for all i. We will define 
two different ensembles: 

1. the ensemble Sx,n of SAWs of length A^ such that ^^(O) = x; 

2. the ensemble of SAWs of length A^ such that both endpoints can be anywhere in 
the cylinder. 

^Note that in a bilocal move a large piece of the walk moves of one or two lattice steps. In protein models 
with nearest-neighbour interactions this may generate large energy differences. In this case, it is probably 
important to consider interactions with a longer range. This corresponds to considering microscopic motions 
of the chain that are small compared to the interaction range. Of course, these problems do not arise in 
continuum models, in which one can make steps as small as one wishes. 

^ In two dimensions one often calls cylinder a strip with periodic boundary conditions. Note that the 
definition given here is different. 



4 



(A) 



(B) 



(C) 



Figure 1: All one-bead moves. (A) One-bead flip. (B) 90° end-bond rotation. (C) 180" 
end-bond rotation. 



(D) 



(E) 



Figure 2: Kink rotations or crankshaft moves. 



In free space the two ensembles are equivalent as long as one is interested in properties of 
the walk itself. In confined geometries they are different since translation invariance is lost 
m. d — D directions. Both ensembles are of physical interest: in a slab one can study the 
statistical properties of polymers that can move freely between the confining surfaces, or 
one can determine the behaviour of polymers grafted at one of the boundaries. 

We wish now to define some algorithms that sample these ensembles of walks. They 
use local and bilocal moves. A local move is one that alters only a few consecutive beads 
of the SAW, leaving the other sites unchanged. A bilocal move is instead one that alters 
two disjoint small groups of consecutive sites of the walk; these two groups may in general 
be very far from each other. 

In our study we will introduce three types of local moves (see Figs. |^, |^): 
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[LO] One-bead flips in which one bead only is moved. 

[LOO] Kink rotations (also called crankshaft moves) in which a three-step kink is rotated. 

[LI] End-bond rotations in which the last step of the walk is rotated. The same move can 
also be applied to the first step of the walk. 

We will then introduce several types of bifocal moves: 

[B22] Kink-transport moves in which a kink is cleaved from the walk and attached at two 
other different sites (see Fig. |^); note that the new kink is allowed to occupy one or 
both of the sites abandoned by the old kink. 

[BKE] Kink-end and end- kink reptation moves (see Fig. ^). In the kink-end reptation move 
a kink is deleted at one location along the walk and two new bonds are appended in 
arbitrary directions at the free endpoint of the walk. Viceversa, an end-kink reptation 
move consists in deleting two bonds from the end of the walk and in inserting a kink, 
in arbitrary orientation, at some location along the walk. The same move can also 
be applied to the first step of the walk. 

[BEE] Reptation move (see Fig. ^) in which one bond is deleted from one end of the walk 
and a new bond is appended in arbitrary direction at the other end. 

We wish now to define algorithms made up with the moves we have presented above 
and that are ergodic. As shown by Madras and Sokal ||^, there exists no ergodic algorithm 
made up of local moves. It is therefore necessary to add some bifocal moves to obtain 
ergodicity. The oldest bilocal algorithm is the reptation algorithm, which uses only the 
moves BEE. As it was soon realized, it is not ergodic due to the presence of walks with 
trapped endpoints. Here we wish to define new bilocal algorithms. We will first consider 
the ensemble £x,n- We will discuss two algorithms: 

• Kink-kink bilocal algorithm. It uses the local moves LO, LI, and the bilocal moves 
B22. The local move LI is applied only to the last point of the walk, otherwise uj{0) 
would not be kept fixed. We will show that it is ergodic in two dimensions (under 
some technical conditions), and that it is not ergodic in three dimensions due to the 
possibility of knots. In higher dimensions its ergodicity is an open problem. 

• Kink- end reptation. It uses the moves BKE applied to the last step of the walk. We 
will show that it is ergodic in a D-dimensional cylinder for d > 3. In two dimensions 
it is ergodic only in free space or in the presence of a single surface, i.e. for D = 1 
and W2 = oo. There exists an extension that is ergodic in a two-dimensional strip: it 
uses the moves BKE applied to the last step of the walk and the local moves LO. 
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It is trivial to modify these algorithms so that the first point is not kept fixed. It is enough 
to apply LI and BKE moves to the first step of the walk, too. However, these modifications 
are not ergodic in the ensemble £n'- indeed, LI and BKE moves never change the parity of 
the first point. If there is translation invariance in one infinite direction, this limitation is 
irrelevant: all walk properties can still be obtained correctly from the ensemble of walks in 
which the parity of is fixed. However, this is not the case in the presence of random 
interactions, since translation invariance is completely lost. In order to sample the ensemble 
£n we introduce a different algorithm: 

• Extended reptation. It uses the moves LO, B22, BEE. We will show that it is ergodic 
in two dimensions. The ergodicity for d > 3 is an open problem. In the absence of a 
definite result, an ergodic extension in d = 3 can be obtained by adding BKE moves. 

The ergodicity properties of these algorithms are proved in the next Section. It should 
be noticed that we have not considered the local moves LOO, which are not necessary for 
ergodicity, but that can be added to the moves LO if one wishes. 

The bilocal moves we have introduced above have already been discussed in the litera- 
ture. Removals and insertions of kinks were introduced in [p| -|43[| in order to study cyclic 
polymers and polymers with fixed endpoints with varying length A^. An ergodic algorithm 
was introduced by Reiter [^: he considers moves LO, LI, B22, BEE, BKE and proves the 
ergodicity of the dynamics (his ergodicity proof requires only BKE and BEE moves) in free 



space. A similar algorithm was used in |^6[ in a study of cyclic polymers, considering LO, 
LOO, and BEE moves. In two dimensions the algorithm is ergodic even in a strip, as we 
shall show below, while in three dimensions it is not ergodic since it does not change the 
knot type of the loop. It is unknown if it is ergodic in a given knot class. 



3 Ergodicity 
3.1 Definitions 

In this Section we will prove various ergodicity theorems for the algorithms we have intro- 
duced before. 



Let us begin by introducing some definitions, following Ref. |f44|| : 

Definition 1: A subwalk ^^[i, j], {0 < i < j < N) is a C-turn of u if j — i > 3, uj[i + l,j — l] 
lies on a straight line that is perpendicular to the steps Auj{i) and Auj{j — 1) and Auj{i) = 
-Auj{j - 1). 

The length of the C-turn is the length of the subwalk. 
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Figure 3: The kink-transport move. A kink has been cleaved from AB and attached at 
CD. Note that the new kink is permitted to occupy one or both of the sites abandoned by 
the old kink. 




Figure 4: The kink-end reptation ( — >) and end- kink reptation {< — ) moves. In ( — >), a 
kink has been cleaved from AB and two new steps have been attached at the end marked 
X. Note that the new end steps are permitted to occupy one or both of the sites abandoned 
by the kink. 
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Figure 5: The reptation move. The head of the walk is indicated by X. The dashed hnes 
indicate the proposed new step and the abandoned old step. 

We say that a C-turn belongs to a line (or surface), or that a line contains a C-turn, if the 
segment a;[i + 1, j — 1] lies on the line (or surface). 

Definition 2: A C-turn of uj, u![i,j] is obstructed if there is a site of uj lying on the line 
segment whose endpoints are a;(i) and u!{j). Otherwise it is unobstructed. 

Definition 3: The enveloping hyper-rectangle R[ijj] of the walk u is the lattice hyper- 
rectangle of minimal volume (which may be zero) containing u. A lattice hyper-rectanglc 
is the set of points {xi, . . . Xd), such that h < Xi < Li, I2 < X2 < L2, . . ., < a;^ < L^, for 
some integers li, I2, ■ ■ ., Li, L2 . . . 

Definition 4' A walk cu is directed if there are no steps that have opposite directions. 

Definition 5: A tower of links of height /i > is a subwalk c<j[i, j] with j — i = 2h + 1, 
< i < j < N, such that u![i, i-\- h] and -|- 1 -|- h,j] are segments and uj{i) and Lj{j) are 
lattice nearest neighbours. 

We call the lattice hnk I connecting u;{i) and u;{j) the base of the tower. We denote the 

tower as T(/, h). 

Moreover we will say that a tower is parallel to a given line if the segments u![i,i + h] and 
-|- 1 -|- h,j] are parallel to this line. 

Definition 6: Given a walk u and a tower T{1, h), I connecting uj{i) and ui{j), we define the 
quotient walk uj/T^l, h) as the walk with sites a;(0) .... a;(i), u!{j), .... u!{N). 

Definition 7: Given a walk uj with a tower T{1, h), let UJ — uj/T{l, h) be the quotient walk. 
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If uj is directed, uj is said to be quotient- directed. 



For the two dimensional proofs we will make extensive use of the following theorem due 
to Madras and reported in Theorem 9.7.2, pag. 356: 

Theorem 1: In two dimensions, if a walk uj has at least one C-turn, then u has an unob- 
structed C-turn. 



3.2 Ergodicity properties of the pivot algorithm in a confined 
domain 

In this Section we want to discuss the ergodicity properties of the pivot algorithm. Following 



Ref. [|T^] it is possible to prove that the algorithm is ergodic in the presence of a single 



confining surface. More precisely, the following theorem holds: 

Theorem 2: In d dimensions, the pivot algorithm is ergodic in a (d — l)-dimensional cylinder 
if < msix{uJd{0),Wd — ujd{0)), where ujd{0) is the d-th component of uj{0) and is the 
number of steps of the walk. 

However, the algorithm is not expected to be ergodic in more constrained geometries. 
We will now show that the pivot algorithm is not ergodic in a two-dimensional strip of 
width w, for N > {w + 1)^. Indeed, pick a bead uj{i) and let d be its distance from the 
boundary y = 0. Then, consider the reflections with respect to the diagonals and the ±90° 
rotations. These are the only transformations that change the number of links that are 
oriented in the ±x and in the ±y directions. It is easy to see that these moves are successful 
only if either —d < uj{j)x — uj{i)x < w — d or d — w < uj{j)x—u!{i)x < d for all j > i. But this 
cannot be verified if N — i > {w Consider now the subwalk Q = uj[0, N — {w + 1)^]. 

The previous argument shows that the number of links belonging to Q that are directed in 
the ±y or in the ±x is fixed. Thus, the algorithm is not ergodic. We believe, although we 
have not been able to prove, that the algorithm is also not ergodic in a three-dimensional 
slab. 



3.3 Ergodicity of the kink-kink bilocal algorithm 

We will now prove the ergodicity of the kink-kink bilocal algorithm in two dimensions in a 
strip of width w = W2- To simplify writing the walks, we indicate by N and S the positive 
and negative y-direction, and by W and E the positive and negative x-direction. Let us 
begin by proving the following lemmas. In all cases we assume d = 2. 
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Lemma 1: Consider a directed walk u and suppose that the distance between cj(0) and at 
least one boundary of the strip is larger than or equal to 2. Then cu can be reduced to any 
given rod using the kink-kink bilocal algorithm. 

Proof : Using LO and LI moves, it is trivial to show that it is possible to reduce the walk 
either to or to . We will now show that it is possible to deform one rod into the 
other, li N <2, the procedure is trivial and thus we will assume N >3. 
Using LO and LI moves, we can deform E^ as follows: 

^jv ^ E^'-^N NE""-' NE'^-^N N^E""'^ (1) 

where we have assumed that the distance between a;(0) and the upper boundary of the 
strip is at least 2. If this is not the case, by hypothesis, the distance between a;(0) and the 
lower boundary of the strip is at least 2, so that the rod E^ can be analogously reduced to 
^2^Ar-2 rjj-^g steps we will present below should then be changed replacing N by S. Then, 
by repeatedly performing the following sequence (p > 2) 

W'^N^EP W^N^EP-'^S W'^N'^ESEP-^ W''+^NEP 

W'^+^NEP-^N W''+^N^EP-\ (2) 

we deform the walk into W^~^N'^E. Finally 

W^'^N^E W^'^N^W W^-'^N^ W^. (3) 

Q.E.D. 

The proof of this lemma requires a technical hypothesis, which however is relevant only 
if w < 2. For larger strips it is always satisfied. In the following we will also assume that 
a;(0) is not a nearest neighbour of a boundary of the strip. Thus, the hypothesis of the 
lemma is satisfied also for w = 2. 

Lemma 2: Consider a quotient-directed walk a; in a two-dimensional strip of width w > 2 
and suppose that a;(0) is not a nearest neighbour of a boundary of the strip. Then uj can 
be reduced to any given rod using the kink-kink bilocal algorithm. 

Proof : Consider uJ — oj/T{l, h), and assume that all steps of u are directed in the A'^, E 
directions. It is immediate to verify that, by using LO and LI moves, one can modify the 
walk obtaining one of these four possibilities: 

(a) N^^E^^N^ES'^E^^- 

(b) n^^e^^s^en'^e''^- 

(c) E'^^N^^E^NW^NPE^^- 

(d) E^^N'^^W^NE^NPE'^K 
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Of course we assume h > 0; otherwise the walk is already directed and it can be re- 
duced to any rod by lemma 1. In case (a), wc can use LO moves to modify the walk into 
j^k^+h-^k2+k3+igh g^^^ ^Yien, combining LO and LI moves, into N'''+^E''+''^+''^+\ The new 
walk is directed and, by lemma 1, it can be reduced to any given rod. 
In case (b), using LO and LI moves, we modify the walk as follows: 

jV^'i E'^^S'^ E^^^^N^ — > N^^ E^'^S^ £;*:3+i jy/i-i ^ _^ jyfci ^fe^/i ^ki+2^h-i _^ 

J\jki j^k2 gh j^h+kz+l ^ ^ j\^ki £<h+k2+ks+l gh ^ jyki £<h+k2+k3+l gh-1 £j ^ ^ 

J\j-ki £<h+k2+k3+2 ^h-1 ^ ^ jyki £<2h+k2+k3+l 

which is directed. By lemma 1, it can be reduced to any given rod. 

Let us now consider case (c). If p > 1, using LO and LI moves, we can deform the walk 
into a new one with p = 1. If p = 1 and > 0, using B22 moves we obtain 

^kij^k2 E^NW^NE^^ E^^N^'^ E^'^'^ NW^ N E^^~^ . . . ^ E^^N^'^ E^^^^NW^N. (5) 

and then, using 

Eki]\[k2 E^^^^NW^N E^^N^"^ E^^^^NW^^^ . . . ^ 

■^ki+ks+h jyl+k2-^h+l ^ ^ki+ks+h jy^h+k2+l 

It is obvious that these last transformations can also be applied when p = (in this case 
we have also — 0). Therefore all walks can be transformed into new ones of the form 
E^NW^. If g = the walk is directed, while for g = 1 wc can transform it into E^NE 
which is also directed. For q > 2 we transform the walk into E^N'^W'^'^ . Since uj{0) is 
not a nearest neighbour of the boundary there is no obstruction to this transformation. If 
g = 2, by means of an LI move, we obtain a directed walk. For g > 3 we apply repeatedly 
the following sequence 

EPN^W"-^ EPN^W-'^S ^ E^N^WSW^-^ ^ 

EP+^NW'^-'^ EP+^N^W^-^, (7) 

obtaining a walk with q — 2 and then a directed walk. Using lemma 1, the walk can be 
reduced to any rod. 

Finally we consider case (d). If k2 = 0, since h > 0, we have ki = 0. Then using LO 
moves we can deform the walk into W^N^'^^E'^'^'^^. If k2 > 0, using LO moves we can modify 
the walk into E^^ NW^ N^'^'^p E^'^^^ . Then, by applying repeatedly the transformation 



E^^~'^NW'^NEN^'^^P~'^E^^^^ E^^''^ NW'^ N^'^^P E^^^'^^'^ , (8) 
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we obtain a new walk of the form NW^ N''^+p E''^^^^+^ and then N''^+p+^ E^^-^^^-^^ . Thus 
the walk can always be deformed into a new one of the form W^N^'E'^, which can be reduced 
to a rod following the method applied to E^NW^. 
Q.E.D. 

It should be noticed that the hypothesis of the lemma is a necessary condition for its 
validity. Indeed, consider a walk with cj(0) = (0, 1), a;(A^ - 2) = (-1, 0), cu(A^ - 1) = (0,0), 
ti;(A^) = (1,0). This walk cannot be deformed into a rod, since the endpoint is frozen. 
This hypothesis is not required if one considers walks for which the first point is not fixed, 
allowing moves LI also on the first step of the walk. 

Lemma 3: Consider a walk in a two-dimensional strip of width w > 2. Then u can be 

reduced to a quotient-directed walk using the kink-kink bilocal algorithm. 

Proof: 

Let us first introduce a few notations. To every walk uj (which is not a rod) with tower 
T(/, h) we associate a triple (tu, T(/, h),t). To define t, set cJ = uj/T(l, h) and consider i?[ciJ]. 
Define Ri and R2 as the two sides of R\uj\ which are perpendicular to the boundary of the 
strip. Let us introduce coordinates so that the x-axis is along the strip. If Ri and R2 have 
equations x = xi and x = X2 respectively, and (xq, |/o), (a; at, Un) are the coordinates of ci;(0) 
and uj{N) respectively, it is not restrictive to assume \xi — xn\ < \xi — xo\ and xi < X2- If 
\xi — xn\ 7^ 0, set t = 0. Otherwise find the smallest x such that the line x = x contains a 
step of UJ or the starting point tu(0) (such a line always exists: in the worst case the walk 
is a rod parallel to the strip and x = X2 = xq). Then set t = x — xi. See Fig. || for an 
example. 

The proof is by induction. The inductive step is the following: given {uj,T{l, h),t) such 
that uJ is not directed, the tower is parallel to the strip, all points of T{1, h) lie on the 
PF-side with respect to the base, and u;(0) does not belong to the line x = Xi + t, then there 
is a sequence of moves such that: a. the new walk also has a tower parallel to the strip, 
lying on the W^-side with respect to the base; b. either t increases or t remains constant 
but /i, the height of the tower, increases. 

To start the induction it is sufficient to assume /i = and choose as / an arbitrary step 
perpendicular to the boundary of the strip (it exists since uj is not a rod). At the end of 
the inductive process we will obtain a new walk, which we will continue to name uj, such 
that either uj is directed or a;(0) belongs to the line x = xi + 1. In the first case, we have 
finished. The second case will be discussed at the end. 

Let us now prove the inductive step. Let us first suppose that x = xi + t contains a 
C-turn uj[k, m]. Then, if the base of the tower does not yet belong to oj[k + 1, m — 1], use 
B22 moves to move the tower so that the base is Aa7(a), k+l<a<m — 2, and the tower 
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Figure 6: The definition of t wlien \xi — xn\ = 0. Tlie dotted line lias equation x = x, h is 
the height of the tower and I its base. 

hes on the VF-side of the base. It is easy to see that, thanks to the definition of t, this 
is always possible. If uJ[k, m] is unobstructed, then by applying LO moves it is possible to 
modify uj[k, m] in such a way that the subwalk of uj connecting cJ(a — 1) and uJ{a + 2) is 
a tower directed to W of height h + 1. In this way either t increases by one or t remains 
constant but h increases by one. 

If uj[k,m\ is obstructed, there exists an unobstructed C-turn a7[i, j] (see theorem 1). If 
the base of the tower does not belong to j], reduce it to a kink, cut it and increase the 
height of the tower by one. If it contains the base of the tower, let us notice that uj[k, m] 
must have length at least 4 as it is obstructed. Therefore there exists a step belonging to 
Lj[k + 1, m — 1] which does not belong to uJ[i, j] . Then move the tower on this hnk (this can 
be done simply by B22 moves if this link is not adjacent to the base of the tower, or by LO 
moves if this is not the case) and at this point reduce cj[i,j] to a kink, and then increase 
the height of the tower by one, using B22 moves. In both cases, thanks to the definition 
of t, the tower can always grow in the VT-direction. In this case t remains constant, but h 
increases by one. 

Let us now suppose that the line x — Xi + 1 does not contain any C-turn. Since, by 
hypothesis, it does not contain a;(0), it can contain only a subwalk a7[i, j] such that either 
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j = N — 2h OT uj\j, N — 2h] is a line perpendicular to Ri, lying on the W-side oi x = Xi+t. 
If the base of the tower does not belong to uj[i, j], with an appropriate choice of y-direction 
we have uj[i, N] = S'^~*"*W^*. Then, using LO and LI moves, it is trivial to modify the walk 
into a new one such that u}[i, N] = VT^"*. Therefore t increases of at least \j — i\, proving 
the inductive step. If the base of the tower belongs to cJ[i,j], with an appropriate choice 
of y-direction Lj[i,N] is of the form S''^W''SE''S''''W* with h > 0. Moreover, if xi is the 
x-coordinate of uj{l), we have xi > x for all I < i. We will now distinguish two cases: (1) 
k2 > 0, (2) k2 = 0. 

If k2 > 0, using LO moves, we can rewrite it as S''^~^''^~''^W'^SE'^SWK Then, using B22 
moves followed by local transformations LO, we have 

sm^sE'^sw^ sm'^+^SE^s sm'^+^s'^E''. (9) 

If /i > 1 we can repeatedly modify the walk as follows: 

gq-\^h+tg2^h _^ s^W^^^S^E^'^N S'^W^^^S^ENE^'^ 

gq-[Yh+t+i g _^ S'^W^'^^'^'^SE^~'^S (lo) 

Then we obtain a walk with the original form and h = 1. But a walk S'^W'^S'^E can be 
modified into S'^W'^ S'^W , and then, by means of LO and LI moves, into a rod 
Therefore in case (1) the subwalk uj[i, N] can be transformed into a rod. 

Let us now consider the case — Since /i > 0, we should have t — Q so that 
oo[i,N\ S^^W^SE^ . If the endpoint of the walk docs not belong to the boundary of 
the strip we can repeat the steps presented for case (1). If the endpoint belongs to the 
boundary and ki > 0, we can use LO moves to modify the walk into S''^~^W'''S^E^, which 
can be transformed into a rod as discussed in case (1). li ki — 0, since w > 2, we can apply 
the following transformations: if h — 1 then 

u;[i, N] = WSE WSW W^S W^; (11) 

if h — 2 then 

u;[i,N] = W^SE^ ^ NWS^E ^ NWS'^W ^ 

NW^S'^ > NW^ >W^; (12) 



if /t > 2 then 



uj[i, N] = W^SE^ NWSW^-^SE^-'^ NW^-^S^E''-^ 
NW^-^S'^E^^'^N NW'^-'^S'^ENE^-'^ 

NW^SE''-^ ^ . . . NW'^S'^E^-'^ > NW'^^-^S'^E. (13) 
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Then we can use the same move discussed in the case h = 2, obtaining W'^^^^. Therefore 
in all cases we deform uj[i, N] into a rod. For the new walk the variable t increases at least 
by 2/i + 1 as required. 

The proof of the inductive step has ended. In this way we have shown that in a finite 
number of steps we obtain a walk uj which is quotient-directed or is such that xq = xi+ t. 
In the latter case, if X2 = Xq, uj is L-shaped and thus directed. Thus, we need only to 
study the case X2 > Xq. Again we will proceed by induction. We assume that the walk 
u has a tower parallel to the strip which lies on the E-side with respect to the base. At 
the beginning notice that R2 must certainly contain a C-turn as it contains none of the 
endpoints. Then move the tower on this C-turn. This is always possible. Then we show 
that as long as X2 > Xq and the walk is not directed we can modify it in such a way that the 
tower increases in height. The argument is exactly identical to the one we have previously 
discussed for the case in which a C-turn exists on the line x = Xi + t. Therefore, in a finite 
number of steps, we obtain a walk with X2 = Xq. As we already discussed, this walk is 
quotient-directed. 

Q.E.D. 

It is now trivial to state the ergodicity theorem for the kink-kink bilocal algorithm which 
is a simple consequence of the lemmas proved above: 

Theorem 3: Consider a walk in a two-dimensional strip of width w >2 and suppose that 
ci;(0) is not a nearest neighbour of a boundary of the strip. Then uj can be reduced to a 
given rod using the kink-kink bilocal algorithm. 

The result we presented above applies only to the two-dimensional case. Indeed, the 
algorithm is not ergodic in three dimensions. For instance, consider the walk {N = 18) 

CO ^ {-yf{x){zf{-xf{-zf{xf{y){zf{-x){-yf. (14) 

By direct enumeration, one can verify that it cannot be reduced to a rod. 

The kink-kink bilocal algorithm can also be used to simulate ring polymers. It is enough 
to exclude the LI moves. Such an algorithm was considered in |^,^. It is clear that this 
algorithm is not ergodic in three dimensions since it does not change the knot type of 
the ring. We will now prove its ergodicity in a two-dimensional strip. We will need the 
following lemma: 

Lemma 4'- In two dimensions, consider a closed walk u, with > 4 steps. It contains an 
unobstructed C-turn uj[i,j] with < i < j < N. 

Proof: Consider w = uj[0, N — 1]. It is not directed, and thus, by theorem 1, it contains an 
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unobstructed C-turn j] with < i < j < N. Since N > A, it is also unobstructed in cu. 
Q.E.D. 

We can now prove the ergodicity of the algorithm. 

Theorem 4'- Consider a lattice point a; in a two-dimensional strip and the ensemble TZx,n oi 
A^-step self-avoiding polygons such that a;(0) = uj{N) = x. The kink-kink bilocal algorithm 
without LI moves is an ergodic algorithm for 71x,n- 

Proof: The proof is similar to that of Lemma 3. Consider a walk cu with tower T{1, h) along 
the strip and let uJ — oj/T{l, h). 

We will now prove the following: if the length of uj is larger than 4, we can modify the 
walk into a new one uo' with tower T(/', h') parallel to the strip and h' > h. Therefore, in a 
finite number of steps, ZJ is reduced to a square of length 4 and a; is a rectangle of height 
one. It is trivial to show that all these rectangle can be modified one into the other ending 
the proof. 

To prove the previous statement, first notice that at least one of the boundaries of the 
enveloping rectangle R[uj] perpendicular to the strip does not contain uj{0) and contains 
a C-turn uJ[k,m], < k < m < N. Using B22 moves we first move the tower on a link 
I e uJ[k + l,m — 1], adding kinks on top of this link, outward with respect to R\uJ]. If 
uj[k, m] is unobstructed we can increase the height of the tower using LO moves, otherwise, 
by lemma 4, there exists an unobstructed C-turn j] such that cj(0) ^ uj[i + — 1]. As 
we already discussed in Theorem 3, such a C-turn can be reduced to a kink which is then 
moved on top of the tower, increasing its height. 

Q.E.D. 

3.4 Ergodicity for the extended reptation algorithm 

We will now prove the ergodicity in two dimensions of the extended reptation algorithm. 
We do not know whether this algorithm is ergodic for larger values of d. 

Theorem 5: The extended reptation algorithm is ergodic in d = 2, for w = W2 > 1. 
Proof: The proof is similar to that of Lemma 3. We assume we have a walk uj with a tower 
T{1, h) parallel to the x-axis. Then we consider the enveloping rectangle i?[aJ]: Ri and i?2 
are the two sides of i?[aJ] perpendicular to the boundary of the strip. We will now show 
that it is possible to deform the walk so that one of the endpoints belongs either to Ri or 
i?2 (defined for the new walk). 

The proof is by induction: given a walk u such that none of the endpoints belongs to 
Ri or i?2 and which has a tower T(/, h) parallel to the x-axis, we can deform it into a new 
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walk with tower T{1', h+ 1) parallel to the x-axis. To prove this statement, notice that Ri 
must contain a C-turn cJ[i,j]. If / ^ uj[i + l,j — 1], using B22 moves, we can modify the 
walk so that it has a tower of height /i on a link belonging to uj[i + 1, j — 1] on the W-side 
with respect to Ri. Since the walk has a C-turn, by theorem 1, it has an unobstructed 
C-turn. If a;[i, j] is unobstructed, we can reduce it to a kink obtaining a tower of height 
h + 1. If uj[k, I] with / < i or A; > J is unobstructed, we can reduce it to a kink and use B22 
moves to move the kink on top of the tower increasing its height. We must then consider 
the special cases in which cJ[i, j] r\uj[k, I] is not empty and cJ[i, j] is obstructed. It is easy to 
realize that one should have l = i + 2 or k = j — 2 and that they can be treated similarly. If 
I — i + 2, since \j — i\ > 4 j] is obstructed) we can modify the walk (if needed) so that 
the base of the tower does not belong to u;[k, I]. Then we reduce it to a kink, and move it 
on top of the tower increasing its height, li k = j — 2 we can proceed analogously. We have 
thus proved the inductive step and therefore we have shown that in a finite number of steps 
we can obtain a new walk such that: (a) one of the endpoints belongs to i?i or R2; (b) 
the base of the tower belongs to Ri. If one of the endpoints belongs to R2, using reptation 
moves in the E-direction, we can reduce the walk to a rod. If none of the endpoints belongs 
to R2, an endpoint belongs to Ri and R2 contains a C-turn cl;[i, j]. Then, by means of B22 
moves, one can move the tower on one of the links belonging to u;[i, j], so that all points of 
the tower lie on the E-side with respect to the base. Then, consider the endpoint belonging 
to Ri. Using reptation moves in the W-direction, we can reduce the walk to a rod. Finally 
note that, ii w > 1, the two rods and can be deformed one into the other. We 
have thus proved the ergodicity of the algorithm. 
Q.E.D. 

It is easy to see that the moves that are added to the reptation algorithm (LO and B22) 
are necessary for the ergodicity of the algorithm. Indeed consider the walk with N — 22: 

N^W^S^ESE^N^E^SESVN^. (15) 

It is easy to see that it cannot be deformed without using LO moves. The walk with N — 13 

NWS^E^NE^S^WN (16) 

requires instead the B22 moves to be reduced to a rod. 

3.5 Ergodicity for the kink-end reptation algorithm 

In this Section we consider the kink-end reptation algorithm. We show that this algorithm 
is ergodic in any cylinder in c? > 3 and in the presence of a single surface (or in free space) 
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Figure 7: With the definitions given in Theorem 6, the dotted vertical fines fiave equations 
X = x^. Tfie step in boldface on tfie left is an S'-step, wfiile tfie step in boldface on tfie 
rigfit is a step belonging to x = x~^ tfiat does not fiave an empty sfiadow and is not an 
^-step. 

in two dimensions. Tfiis algoritfim is not ergodic in a two-dimensional strip. However 
ergodicity is recovered by adding LO moves. 

Theorem 6: For c/ > 3, tfie kink-end reptation algoritfim is ergodic for w = miUjWj > 1 
and > 3. 

Proof: In order to prove tfie tfieorem let us introduce some useful definitions. We say tfiat 
a walk step Au^j), 1 < j < — 4, is an S-step (see Fig. |^) if Auj{j — 1) = Au^j + 1) 
are directed in tfie (±a;)-direction, wfiile Auj{j) is ortfiogonal to tfiem. Note tfiat fiere we 
indicate by x tfie first direction tfiat, by definition, is always infinite. Given a walk step 
Auj{j) wfiicfi is not in tfie x-direction, we also define its positive and negative shadow. If 
X is tfie unit vector (1, 0, . . . , 0), tfie positive (resp. negative) sfiadow is tfie set of lattice 
points {uj{j) + nx,uj{j + 1) + nx}, wfiere n is a positive (resp. negative) integer. We say 
tfiat Auj{j) fias an empty shadow if its positive or negative sfiadow contains none of tfie 
walk sites a;(/c), < k < N — 2. To understand tfie relevance of tfiis definition, note tfiat, if 
ALo'(j) fias an empty sfiadow, tfien (see Fig. |^) we can perform end- kink moves by adding 
kinks on top of Auj{j), obtaining eventually a new walk uj' sucfi tfiat — 2] is a rod in 
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the x-direction. Note that, by definition, the last two walk sites may belong to the shadow. 
However, they do not represent an obstruction for the end-kink moves since the last two 
steps are deleted in the first iteration of the process. 

The proof of the theorem is based on the following inductive step: given a walk uj such 
that the subwalk uj[i, A^ — 2],0<z<A^ — 2, isa rod in the (±x)-direction, we can deform 
it into a new walk uj' such that — 2], j < z, is also a rod in the (±x)-direction. Note 

that it is possible that, at the beginning of the induction process, i = N — 2: it simply 
means that Auj{N — 3) is not in the (±x)-direction. This inductive step allows to prove 
that, in a finite number of steps, any walk can be deformed into {±x)^~'^X where X is the 
subwalk uj[N — 2, A^]. Then, since w > 1 and N > 3, with an appropriate choice of the 
y-axis, we can deform it as follows: 

{±xf''X y{±x)i-y)i±xr-' ^ (±x)^. (17) 

Thus all walks can be deformed into a rod (±x)^. Finally, it is easy to show that the rod 
can be deformed into (— x)^, proving the ergodicity of the algorithm. 

To prove the inductive step, let us introduce coordinates u{k) = {xk,yk, Zk, ■ ■ ■) and 
define x~ (resp. 5;+) as the smallest (resp. largest) value of x such that there exists a walk 
site uj{k), k < N — 2, with Xk = x'^ and Auj{k) not in the (±a;)-direction. If uj[0, — 2] is 
not a rod in the (±a;)-direction, x~ and x~^ certainly exist although they may coincide. 

Now consider the links belonging to the hyper-surfaces x = x~ and x = W^. Suppose 
that one of them Aa;(/), / < A^ — 2, is not an S'-step. It is not restrictive to assume that it 
belongs to x = x^. We will now show that there is Auj{j), contained in x = x~^ that has 
an empty shadow, with j < i. If Aci;(/) has an empty shadow, we can take j = I. Since 
I < N — 2, and all steps Auj{i), . . . Auj{N — 3) are in the x-direction, we have j < i. If 
Aa;(/) does not have an empty shadow, there are two possibilities: (a) Auj{l — 1) is in the 
negative x-direction; (b) Auo{l + l) is in the positive x-direction. We will consider only case 
(a) since case (b) is completely analogous. In case (a) (see Fig. |^ we will now show that 
Au}{l + 1) has an empty shadow and that / + 1 < i, so that we can take j = l + l. To prove 
this statement we will show the following: (al) / < A^ — 3; (a2) Ati;(Z + l) belongs to x = x"^; 
(a3) Au}{l + 2) cannot be oriented in the positive x-direction. /^From (a2) and (a3) we see 
that Aci;(/-t-l) has an empty positive shadow, while (a2) and (al) allow to conclude l + l < i 
as required. To prove (al), note that, if (al) were not true, we would have I = N — 3 and 
the walk would be of the form {—x)^~^diX , where X = uj[N — 2, N]. But this implies that 
Auj{l) has an empty (negative) shadow, which is against the initial assumption. To prove 
(a2), note that Auj{l + 1) cannot be oriented in the negative x-direction; otherwise, since 
/ < A^ — 3, Au;(/) would be an 5'-step. If it were directed in the positive x-direction, then 
the walk would be (— x)'c?x^~'~^X where d is the direction of Auj{l) and X indicates the 
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last two steps. But in this case Au{l) would have an empty (negative) shadow. Therefore, 
(a2) is proved. If (a3) were not true, the walk would be (— x)'(ii(i;+iX^~'~'^X, where di and 
di+i are the directions of Aa;(/) and Auj{l + 1) respectively, and, therefore, Alj{1) would 
have an empty (negative) shadow, against the initial hypothesis. 

If Au{j) has an empty shadow, then we can deform the last steps of the walk as follows 
(we assume, without loss of generality, to have a positive empty shadow): 

uj[j, N] xy{~x) . . . ^ x^y{-xy ... ^ ... ^ x^-^-^y{-x), (18) 

where y is the direction of Au{j) (see Fig. ||). We have thus proved the inductive step. 

Let us now suppose that all walk steps belonging to a; = x" and x = x^ are S'-steps. 
In this case it is easy to convince oneself that, with an appropriate choice of axes, the walk 
has the form x^yXx^Y , where j > 0, /i > 0, X is the subwalk uj[j + 1, N — h — 2] and Y 
is the configuration of the last two steps. Clearly, if uj{k) G X, then < < x"*". If 
x~ = x"*", X is empty and uj = x^yx^Y. Consider now Auj{j) which is the only walk step0 
(as it can be seen from the explicit expressions above) belonging to the face x = x". Since 
w > 1, it is possible to fix the positive z-direction in such a way that uj{j) + z {z is the 
unit vector in the positive 2;-direction) is inside the cylinder. Then we deform the walk as 
follows: 

UJ zx{—z)x^~^ . . . — * {—x)zx'^{—z)x^~^ . . . ^ (— x)^~^2;x, (19) 

which can be reduced to a rod. Thus in all cases we have proved the inductive step. 
Q.E.D. 

It is easy to see that this algorithm is not ergodic in a two-dimensional strip of width 
w = W2- Consider for instance the walks of the following form: 

{{ENYE{ESYEfE'^ (20) 

with k>l. The length of such a walk is X = (4w + 2)A; + 2. It is easy to verify that they 
cannot be modified by the algorithm. 

The previous theorem can be extended to two dimensions in free space, or in the presence 
of a single boundary, i.e. for w = oo. 

Theorem 7: In two dimensions the kink-end reptation algorithm is ergodic for X > 3, in 
the presence of a single boundary or in free space. 

^Note however that <jj{N) may belong to the face x — x~ . This is of no relevance for the following 
discussion since in an end-kink move the last two steps are deleted. 
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Proof: The proof follows the previous one. We should only modify the proof of the 
inductive step in the case in which only S'-steps belong to the lines x = x~ and x = x^ . In 
this case, with a proper choice of axes, the walk has the form 

uj = EPNXE^Y, (21) 

p > 1, h > 1, where X is the subwalk uj[p + 1, N — h — 2] which is contained between the 
lines X = x^ and Y = uj[N — 2, N]. If uj{p) does not belong to the boundary, we can modify 
the walk as follows: 



UJ 



E'PNXE^Y SENE'P-^NXE^ ^ WSE'^NE^'^ ...^ W^-'^SE, (22) 



which can be reduced to a rod. 

We should finally discuss the caseQ in which ct;[0,p] belongs to the boundary y = 0. Let 
y be the largest value of y such that there exists a walk step Acij(Z), I < N — 2, belonging to 
the line y = y- First we show that, iiyk>y {yu is the y-coordinate of uj{k)), then k = N—1 
OT k = N. Indeed, since yo = 0, if we had yk > V and k < N — 2, then u![k — 1, — 2] 
would point N which is in contrast with (^T|). Now, consider the sites belonging to the 
line y = y and let uj{k) be the site with smallest k. Clearly k < i. Note, moreover, that 
Auj{k — 1) is in the positive y direction. Indeed, it does not lie on y = y, otherwise uj{k — 1) 
would lie on this line (remember that uj{k) is the walk site with smallest k belonging to 
this line). It is not oriented in the negative y-direction, otherwise, yk-i > y- For the same 
reason Auj{k + 1) is not in the positive y-direction, unless k + 1 = N — 2. Then, we can 
use end- kink moves to put kinks on top of Auj{k) in the positive y direction. There is no 
obstruction to these moves, since the only possible walk sites that can have a larger y are 
the last two sites that are removed in the first iteration of the process. In this way we 
modify the walk into a new one (we keep calling it u) such that uj[k — 1, — 2] is a rod 
directed in the positive y-direction. 

Now let X be the largest x such that there exists a walk step Ati;(/), / < N — 2, belonging 
to the line x = x. Consider the steps belonging to this line and let j be the smallest integer 
such that Auj^j) belongs to x = x. Clearly j<k — l<k<i. Since ci;[0,p], p > is 
directed E and — 2] is directed N, Auj{j) has an empty positive shadow. Therefore, 

by means of end-kink moves we modify the walk so that — 2] is a rod in the positive 

direction. Since j < i, we have proved the inductive step. 

Q.E.D. 

If one considers two dimensional strips, an ergodic algorithm can be obtained adding 
LO moves. 



^li p > 1, we can use the sequence of moves (22) adding kinks "above" the boundary. The difficult case 
corresponds to p = 1. 
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Theorem 8: In two dimensions the kink-end reptation algorithm with LO moves is ergodic 
for > 3 and W2 > 1- 

Proof: The proof is identical to that of Theorem 6. We should only change the proof of 
the inductive step for the case in which there are only S-steps on the lines a; = x^. With 
a proper choice of axes the walk has the form (0). Then, using LO moves, we modify the 
walk as follows 

UJ = EP-^NE ... ^ ... ^ NEP .... (23) 
Then, using end-kink moves, 

u WNEP^'^ ... ^ ... ^ W^-^NE, (24) 

which can be deformed into a rod. 
Q.E.D. 

Although the kink-end reptation algorithm is ergodic in two dimensions in the absence 
of confining surfaces, the proofs of the theorems indicate that "staircase" sections of the 
walk (for instance sections of the form . . . ENENENEN . . .) will be changed very slowly 
by the algorithm. Therefore, in order to have an efficient implementation, it is probably 
useful to include in all cases the LO moves. 

4 Transition matrices 

In the previous Section we have discussed the ergodicity of the algorithms. Now, we discuss 
how to implement them in order to obtain the correct probability distribution. Here we 
will discuss how to use them to generate walks with uniform probability in the ensembles 
Sx^N or S^. Any other probability distribution can be obtained by adding a Metropolis test 
or a generalization thereolQ. 

4.1 Kink-kink bilocal algorithm 

We will begin by considering the kink-kink bilocal algorithm. Although not necessary to 
ensure the ergodicity of the algorithm, we will also add the local LOO (crankshaft) moves. 
In order to describe the algorithm it is important to classify the possible configurations of 
three successive links (see Fig. |^): 

1. the bonds have the same direction (I configuration); 

^ For simple lattice models of homopolymers and proteins, the Metropolis criterion should provide a 
plausible physical dynamics |45j. For a different point of view see ^ , |47|] . 
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2. two consecutive bonds have the same direction, while the third one is perpendicular 
to them (L config.); 

3. the first and the third bond are perpendicular to the second one, and they are either 
parallel or perpendicular to each other (S config.); 

4. the first and the third bond are perpendicular to the second one, and they are an- 
tiparallel to each other (U config.). 

The algorithm works as follows: 

• Step 1. Choose a random site i of the current walk cu, < i < N. li i — N, propose 
an LI move and go to step 5. 

• Step 2. Determine the configuration of the subwalk uj[i — l,i + 2]. If i = N — 1 we 
imagine adding a link Auj{N) parallel to Au!{N—l), so that the possible configurations 
are of type L and I. Analogously if i = 0, we imagine adding a link Aa;(— 1) parallel 
to Acc;(0). 

• Step 3. Depending on the configuration of a;[i — 1, i + 2], do the following: 

1. I: with probability 1 — {2d — 2)p(22) we perform a null transition. Otherwise we 
go the next step. 

2. L: with probability 1 — {2d — 3)p(22) — p(0) we perform a null transition, with 
probability p{0) we propose an LO move and go to step 5. Otherwise we go to 

the next step. 

3. S: with probability 1 — {2d — 4)p(22) — 2p{0) we perform a null transition, with 
probability 2p{0) we try an LO move: in this case there are two possibilities and 
we choose among them with equal probability and then go to step 5. Otherwise 
we go to the next step. 

4. U: with probability {2d — 3)p(00) we try an LOO move: there are {2d — 3) pos- 
sibilities which are chosen randomly; then we go to step 5. Otherwise we go to 
the next step. 

• Step 4. Choose a second integer j uniformly in the disjoint intervals, — 1 < j < N, 
j ^ i — + 1. If J = — 1,A^ make a null transition. Then, depending on the 
configuration of — 1, i + 2], do the following: 
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— u![i — l,i + 2] is of type I, S, L: if j' = or j = A?" — 1, or if uj[j — l,j + 2] is not 
of type U perform a null transition. Otherwise propose a B22 move, cutting the 
kink u[j — l,j + 2] and adding it to a;[i, i + 1] in one of the possible directions. 
Then go to the next step. 

— a;[i — 1, i + 2] is of type U: according to the configuration of a;[j — 1, j + 2] (if 
j — 0, N — 1 imagine adding links as before) do the following: 

1. a;|j — 1, j + 2] is of type I: with probability (2d— 2)p(22) perform a B22 move: 
cut the kink a; [i — 1 , i + 2] and add it on top of cu [j, j ' + 1] in a possible random 
direction, and then go to step 5. Otherwise perform a mill transition. 

2. uj[j — l,j + 2] is of type L: with probability (2(i — 3)p(22) perform a B22 move: 
cut the kink cc;[i — 1, ^ + 2] and add it on top of a;[j, j + 1] in a possible random 
direction, and then go to step 5. Otherwise perform a null transition. 

3. uj\j — l,j + 2] is of type S: with probability (2d— 4)p(22) perform a B22 move: 
cut the kink a;[i — 1, i + 2] and add it on top oiculj, j + in a possible random 
direction, and then go to step 5. Otherwise perform a null transition. 

4. u}[j — l,j+2] is of type U: with probability {2d—3)p{22) perform a B22 move: 
cut the kink uj[i — l,i + 2] and add it on top of c<j[j, j + in a possible random 
direction, and then go to step 5; with probability {2d — 3)p(22) perform a 
B22 move: cut the kink u[j — l,j + 2] and add it on top of u;[i, i + 1] in a 
possible random direction, and then go to step 5. Otherwise perform a null 
transition. 

• Step 5.: Check for self-avoidance. If the proposed new walk is self-avoiding keep it, 
otherwise perform a null transition. 

The algorithm we have presented depends on three probabilities p(0), p(00) and p(22) 
that are the probabilities of an LO, LOO and B22 move respectively. It is easy to check that 
the algorithm satisfies detailed balance and thus the walks are generated with the correct 
probability distribution. We should now determine the single probabilities that must be 
such to satisfy the obvious constraint 



(25) 



Considering the configurations I, L, and S we obtain the constraints 



{2d-2)p{22) < 1, 
{2d - 3)p{22) + p{0) < 1, 
(2d - 4)p(22) + 2p(0) < 1. 



(26) 
(27) 
(28) 
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li u![i — l,i + 2] is of type U, we obtain, depending on the configuration oi u!\j — 1, j + 1]: 

{2d - 2)p{22) + {2d - 3)p{00) < 1, (29) 

{2d - 3)p{22) + {2d - 3)p{00) < 1, (30) 

(2rf - 4)p(22) + (2rf - 3)p(00) < 1, (31) 

2(2(i-3)p(22) + (2(i-3)p(00) < 1. (32) 

These conditions impose for d > 2: 

{4:d - 6)p{22) + {2d - 3)p{00) < 1, (33) 
(2rf-4)p(22) + 2p(0) < 1. (34) 

A solution of these inequalities which maximizes p(0) and p(00) at p(22) fixed, is 

p(0) = i[l-(2d-4)p(22)], (35) 

= ^^[l-(4d-6)p(22)], (36) 

P(22) < (37) 

Since the LOO move is not necessary for the ergodicity of the algorithm, while the B22 one 
is essential to ensure a fast dynamics, it is natural to require p(22) to be maximal, even if 
this imphes p(00) = 0. Then we obtain the following transition probabihties: 

p{00) = 0, (39) 
P(22) = jj^. (40) 

In two dimensions p(0) = p(22) — 1/2, while in three dimensions p(0) = 1/3 and p(22) = 
1/6. 

4.2 Extended reptation algorithm 

This algorithm extends the standard reptation method. The reptation (or shthering-snake) 
algorithm has two different implementations. The first one, which satisfies detailed balance, 
works as follows: 
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Step 1. With probability 1/2 delete u![N — 1, N] and add a new link at the beginning; 
otherwise delete uj[0, 1] and add a new link at the end of the walk. 



• Step 2. Check if the new walk is self-avoiding. If it is keep it, otherwise perform a 
null transition. 

A second version uses an additional flag which specifies which of a;(0) and uj{N) is the 
"active" endpoint. It works as follows: 

• Step 1. Delete one bond at the "active" endpoint and append a new one at the 
opposite end of the walk. 

• Step 2. If the new walk is self-avoiding keep it, otherwise stay with the old walk, and 
change the flag, switching the active endpoint. 

This algorithm does not satisfy detailed balance, but it satisfies the stationarity condition 
generating the correct probability distribution. 

The extended reptation algorithm consists in combining with non-zero probability the 
reptation algorithm and the kink-kink bilocal algorithm. More precisely the algorithm 
works as follows: 

• Step 1. With probability p perform a reptation move, with probability 1 — p a kink- 
kink bilocal move, as specified in the previous section. 

Note that in this algorithm the LI moves are no longer needed. Therefore one can modify 
Step 1. of the kink-kink bilocal algorithm choosing i such that < i < — 1. The 
probability p is not fixed. It is only required that < p < 1 to ensure the ergodicity of the 
algorithm. It can therefore be tuned in order to obtain the best critical behaviour. 

4.3 Kink-end reptation algorithm 

The kink-end reptation algorithm uses kink-end and end-kink reptation moves (see Fig. 
We will present here two different implementations of the algorithm which, however, are 
expected to have the same critical behaviour. 

Let us explain the first implementation. An iteration consists of the following steps: 

• Step 1. Choose a random site i of the current walk with < i < A^ — 3. 

• Step 2. Try an end-kink move with probability p{EK) or a kink-end move with 
probability p{KE) = 1 ~ p{EK). In the first case delete the last two bonds of the 
walk and insert a kink on the bond Au!{i) in one of the {2d — 2) possible orientations. 
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In the second case, if i 7^ and — 1, i + 2] is a kink, remove it and attach two 
bonds at the end of the walk in one of the {2d — 1)^ possible ways; otherwise perform 
a null transition. 

• Step 3. Check if the proposed walk is self- avoiding. If it is keep it, otherwise make a 
null transition. 

The transition matrix is given by 

if^^EK (41) 

and 

^(^^^')= (]^^P(r)(^^^X5ATy(a;') ifT = KE (42) 
for cu 7^ cu', where xsaw is defined by 

, I 1 if a;' is self-avoiding 
XsAwiu.') - [ (43) 

Detailed balance requires 

p{KE) {2d - 2) = p{EK) {2d -if, (44) 
from which, using p{EK) + p{KE) = 1, we obtain 

pyiyj^) — (2d_i)2+(2d-2)' 

p{EK) = (2d_i)2+(2d_2) > 

In two and three dimensions we obtain explicitly 



(45) 



p{EK) = 2/111 p{EK) = 4/29 1 

p{KE) = 9/11/''''*"^' p{KE) = 25/29 J''''^"'^- ^^"^ 

The second implementation of the algorithm is similar to that of the kink-kink bifocal 
algorithm. An iteration consists of the following steps: 

• Step 1. Choose a random site of the current walk with < i < A'" — 3. 

• Step 2. Determine the configuration of the subwalk — 1, i-|-2]. If i = 0, we imagine 
adding a link Aa;(— 1) parallel to Aa;(0). 
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• Step 3. Depending on the configuration of cufi — 1, i + 2], do tlie following: 

1. I: with probability (2d — 2)p{EK) perform an end-kink move, deleting the last 
two steps and adding a kink in one of the possible {2d — 2) directions; otherwise 
perform a null transition. 

2. L: with probability (2d — 3)p{EK) perform an end-kink move, deleting the last 
two steps and adding a kink in one of the possible {2d — 3) directions; otherwise 
perform a null transition. 

3. S: with probability {2d — 4:)p{EK) perform an end- kink move, deleting the last 
two steps and adding a kink in one of the possible {2d — 4) directions; otherwise 
perform a null transition. 

4. U: with probability {2d — 3)p{EK) perform an end-kink move, deleting the last 
two steps and adding a kink in one of the possible {2d — 3) directions; with 
probability {2d — l)'^p{KE) perform a kink-end move, cutting the kink and 
adding randomly two links to the walk in random directions; otherwise perform 
a null transition. 

• Check whether the proposed new walk is self-avoiding. If it is keep it, otherwise make 
a null transition. 

The transition matrix is given by 

P{u u') = ^^\^^ p{T)xsAw{oo') where T = EK or T = KE (47) 

for Lo 7^ uj', where XsAw is defined in Eg. (^3]) . Detailed balance requires p{KE) = p{EK) = 
p, while Eq. (|25|) gives the following constraints: 

{2d-2)p{EK) <1, (48) 

{2d - 3)p{EK) < 1, (49) 

{2d - 4:)p{EK) < 1, (50) 

{2d - lfp{KE) + {2d - 3)p{EK) < 1. (51) 



It follows 



1 . ^ 

p < — , 52 

^ - 2d-2' ^ ' 

1 

^ - (2rf-3) + (2rf-l)2- 
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For d > 1, the second inequality is the most restrictive one. Therefore the best choice of 
p corresponds to taking the equahty in Eq. (0). In two and three dimensions we obtain 
p = 1/10 {d = 2) and p = 1/28 {d = 3). 



5 Dynamic critical behaviour 

In order to understand the efficiency of an algorithm one should analyze the autocorrelation 
time r. There are several different definitions^ for r: the exponential autocorrelation time 
Texp that controls the relaxation of the slowest mode in the system and the integrated 
autocorrelation time Tint,o that depends on the observable O one is considering and that 
controls the statistical errors on O. For N — »• oo, one expects a dynamic critical behaviour, 
i.e. T ~ N^, where the exponent z may depend on which autocorrelation time one is 
considering. 

We now derive lower bounds on the exponent z in the absence of interactions. Let us 
consider global observables, like the squared end-to-end distance Rl and the squared radius 
of gyration R^. For bilocal algorithms we expect |10|,^ Tmi.o ~ The basic assumption 



is that the slowest mode appearing in global observables is associated to the relaxation of the 
squared radius of gyration Rg. Then, an estimate of r can be obtained as follows. At each 
elementary step, R^ changes by a quantity of order N'^'^~^ . An independent configuration 
is reached when the observable changes by one standard deviation N'^^ . Assuming that the 
observable performs a random walk, we obtain r ~ (A^2i/^jy2i/-i^2 ^ j\j2 ^ practice the 
argument should provide only a lower boundQ which we expect to be correct for all global 
observables. 

In particular, this should apply to the end-to-end distance R^. Here, however, we 
should notice that our algorithms update the end-point of the walk with very different 
frequencies. The extended reptation and the kink-end reptation change ijj{N) every 0(1) 
iterations, while the kink-kink bilocal algorithm updates uj{N) only every 0{1/N) iter- 
ations. Therefore, an additional factor of should be added for the kink-kink bilocal 
algorithm: in this case, we expect r > A^'^. It is interesting to notice that the kink-kink 



^We refer the reader interested in more precise and rigorous statements to [ 44 , 48| . 
In ^] the foUowing rigorous lower bound was obtained: 

2var(C') 1 

where var(C') is the static variance of the observable O and C is the maximum change of O in a single 
Monte Carlo step. In the heuristic argument given above, we have replaced C by the average change of O 
in a single Monte Carlo step. 
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d 


p{\) 


p{l) 


J9(U) 


P(S) 


2 


0.152 


0.481 


0.108 


0.259 


3 


0.051 


0.356 


0.102 


0.491 



Table 1: Probabilities of the different configurations of three links. 



bilocal algorithm behaves approximately as the algorithm of Reiter [^]: indeed, also in 
this case, the endpoint is updated with frequency 1/A^. The available numerical results for 
very short walks (A^ < 100) are in agreement with the bound given above: they indicate 
r ~ A^^ in two and three dimensions. 

While Reiter's algorithm can be easily speeded up by increasing the frequency of the 
BEE and BKE moves, no improvement is possible for the kink-kink bilocal algorithm. 
Indeed, in this case it makes no sense to increase the frequency of the LI moves. Clearly, 
endpoint moves should be performed with the same frequency of the moves that change 
the site u^N — 1), otherwise they do not effectively change the endpoint position. But 
a;(A^ — 1) is updated by B22 and LO moves with frequency 1/A^. Therefore LI moves 
should be performed with the same frequency. 

Let us now discuss in more detail the different implementations of the extended reptation 
and of the kink-end reptation algorithm. For the extended reptation, we should choose 
between the two different implementations of the reptation dynamics. If one considers 
ordinary random walks, it is obvious that a new walk is generated in 0{N'^) iterations of 
the first algorithm and in exactly N iterations of the second one. Thus the second version 
is much more efficient than the first one. For SAWs we do not expect such a big difference 
since the walk will move in a given direction only for a small number of steps 8 in 
two dimensions, ~ 14 in three dimensions in the absence of interactions). Therefore we 
expect an improvement by a constant factor, and, indeed, simulations show that in 
three dimensions the second implementation is 5-6 times faster that the first one. In the 
extended reptation we should also fix the parameter p. From the discussion given above, 
it is clear that we must have p > as A^ ^ oo, otherwise the motion of the endpoints 
slows down the dynamics. To fix its specific value, we may compare p to the probabilities 
of proposing local and bilocal moves in the kink-kink bilocal algorithm. Assuming that the 
probability of occuring of a I, L, U, and S configurations is independent of the position of 
the walk site — it should be approximately true for large values of A^ — the probability pb 
of a bilocal move B22 and the probability pi of a local move LO are given by 

Pb = 2p(22)p(U)[(2rf-2)p(l) + (2t;-3)(p(L)+j9(U)) + (2rf-4)j9(S)], (55) 
pi = p(0)(p(L) + 2p(S)), (56) 
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To have a quantitative prediction we should know p(U), p(L), and p{S). If we were 
considering non-reversal random walk we would have 

P(U) ^ (5T) 



P(L) ^ (59) 
pfS) = ^'''-I'P''-^' (60) 



SO that 



16(d-lY , , 

Pb = rn. o^.o. (61) 



(2d-3)(2rf- 1; 



4(t/- 1)3 

^' ~ (2rf-3)(2d-l)2' ^ ^ 

In two and three dimensions we obtain pb ~ 0.1975, pi ~ 0.4444, and pb ~ 0.1365, ~ 
0.4267 respectively. For SAWs the probabilities can be computed by means of a short 
Monte Carlo simulation. The results are reported in Table |T]. In two dimensions we obtain 
Pb = 0.096, pi = 0.500, while in three dimensions pb = 0.087, pi = 0.446. Thus, in the 
extended reptation algorithm, B22 moves are proposed with probability ~ 0.1(1 —p), while 
reptation moves are proposed with probability p. If one wants to balance these two types of 
moves — this is reasonable of one wants a physical kinetics — one should choose p ~ 0.1. 
On the other hand, it is clear that reptation moves are more relevant than B22 moves. 
Indeed, reptation moves are essential for the motion of the endpoint, while B22 moves are 
required only to avoid the trapping of the endpoints. Thus we expect the algorithm to be 
more efficient for larger values of p. The simulations |^0|] indicate that the fastest dynamics 
is obtained for 0.5 ^ p ^ 0.9. 

Finally let us consider the kink-end reptation algorithm. We presented two different 
versions and we discuss now their relative efficiency. The first implementation chooses 
the move without checking the nearby bonds. A deformation is always proposed but it 
may immediately fail because it does not respect self-avoidance when one considers the 
neighbours of the chosen bond. The second algorithm is more careful: the move is chosen 
after considering the position of the nearby bonds. However, with a finite probability, 
it performs a null transition. In order to compare correctly the two implementations we 
should therefore compute: 
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(a) for the first algorithm, the probabihty of bilocal moves that do not fail after checking 
the position of the two nearby bonds; 

(b) for the second algorithm, the probability of proposing a bilocal move. 

For the first algorithm an easy computation gives 

Pkink-end = p{\})p{KE) (63) 

Pend-kink = [(2d - 2)p(l) + (2d - 3) (p( L) + p(U)) + (2d - 4)p(S)] , (64) 

while for the second we have 

Pkink~end = p(U)p(2d - 1)^ (65) 
Pend-k^nk = P [(2d - 2)p(l) + (2d - 3) (p(L) + p(U)) + (2d - 4)p(S)] . (66) 

Therefore the first algorithm is equivalent to the second one with probability 

_ pjKE) _ p{EK) _ 1 
^ (2d -1)2 (2d -2) (2d -2) + (2d -1)2' ^ ^ 

i.e. p = 1/11, 1/29 in two and three dimensions respectively. This should be compared 
with the optimal value of p, p — 1/10 and p = 1/28 for d = 2, 3. Thus the second 
algorithm is more efficient than the first one, as it should be expected since the second one 
chooses the proposed move more carefully. However the improvement in efficiency is small, 
approximately 10% in two dimensions and only 3% in three dimensions. 
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□ 




Figure 8: The sequence of moves that generate a new walk uj' such that a;'[j, — 2] is a 
rod in the (±a;) direction. The dashed hne in the upper part of the figure indicates an 
arbitrary subwalk. The dotted hne is the hne of equation x — x'^ . 
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(d) 



Figure 9: Configurations of tliree consecutive links: (a) configuration of type I; (b) con- 
figuration of type L; (c) configuration of type S; (d) configuration of type U. 
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